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A large approximate degeneracy appears in the light meson spectrum measured at LEAR, sug- 
gesting a novel principal quantum number n + j in QCD spectra. We recently showed that the 
large degeneracy could not be understood with state of the art confining and chiral invariant quark 
models, derived in a truncated Coulomb gauge. To search for a solution to this problem, here we 
add the gluon or string degrees of freedom. Although independently the quarks or the gluons would 
lead to a 2?! + j or 2n + I spectrum, adding them may lead to the desired n + j pattern. 



To understand the laxge degeneracy of the hght meson 
spectrum [H [3, d, 0, U observed by the Crystal Bar- 
rel collaboration at LEAR in CERN, we include at the 
same token both bosonic and fcrmionic excitations. This 
is relevant for the observation of excited mesons in the 
new generation of experiments, say at PANDA at FAIR, 
and in Lattice QCD. We first describe the degeneracy, 
where the meson mass squared in the excited me- 
son spectrum is led by the principal quantum number 
n + j. Then we review the status of gluon string exci- 
tations with a AP possibly of the order of 2TT<T(2n + I). 
Recently we have shown 5] that the quark degrees of free- 
dom are not sufficient to understand the excited light me- 
son spectrum, finding that the chiral invariant and con- 
fining quark model has a Af^ of the order of 2TTa{2n + j), 
with a pre-factor is correct for the angular / or j quantum 
numbers, but too large for the radial quantum number n. 
Importantly, both gluon and quark excitations have the 
same scale, with the string tension a. Here we show that 
arriving at the desired n + j pattern is possible if we in- 
clude the gluon or string degrees of freedom in relativistic 
quark models with chiral symmetry breaking. For sim- 
plicity, in this first comprehensive study of gluonic and 
chiral invariant high spectrum, we use einbeins to add 
the bosonic excitations to the fermionic excitations. 

In the report of Bugg 0, 0, 13, 0, @ a large degeneracy 
emerges from the spectrum of the angularly and radially 
excited resonances produced in pp annihilation by the 
Crystal Ball collaboration at LEAR in CERN This 
degeneracy is a larger degeneracy than the chiral degen- 
eracy systematically searched by Glozman et al 0, 0, [1] , 
JafFe et al Q, Afonin [l3|; and retrospectively already 
present in the light-light spectrum of Le Yaouanc et al. 

13, [3. 
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FIG. 1: Gluon (left) and string (right) excitations in a meson. 



where ao is the intercept of the il/^ = axis, a is the 
slope and n is a radial excitation. The slope for angular 
and for radial quantum excitations respectively of 
0.877 GeV"^ and 0.855 GeV"^ are almost identical high 
in the spectrum and similar to {2t:(7)~^ — 0.84 GeV~^ 
where the string tension is cr = 0.19 GeV^. This is 
quite remarkable, since in nature there are few examples 
of principal quantum numbers, except for the coulomb 
non-rclativistic potential, say in the hydrogen atom, or 
the harmonic oscillator non-relativistic potential, say in 
the vibrational spectrum of molecules. 

We now discuss the gluon or string degrees of freedom 
in order to include them in the chiral invariant and con- 
fining Quark Model. This is still an open problem, but 
nevertheless in the end of the letter we will see that it 
may lead to the observed degeneracy. The hybrid excita- 
tions of mcsoirs can be addressed with either the bosonic 
Nambu-Gotto-Polyakov string fluctuations or with a con- 
stituent gluon linked to the quarks by fundamental rigid 
strings [16|, as depicted in Fig. [T] Baker and Steinkc [17j . 
only consider two transverse modes of the string with ro- 
tating quarks fixed to the ends of the string, resulting 
in the mass spectrum. 



111 and the heavy-light spectrum of PB et al. 
Also notice that a long time ago, Chew and Frautschi re- 
marked the existence of linear Regge trajectories [ij] for 
angularly excited mesons. A similar linear aligning of ex- 
cited resonances was also reported for radial excitations 
(lit . Presently the status of the excited meson spectrum 
approximately generalizes the Regge trajectories for the 
total angular momentum j. 
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J -t- n ~ ao + aM 
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led by 71 -|- ^ in units of 271 a, similar to the spectrum 
of an open string. This result motivated Afonin Q to 
propose it as the solution to the large n + j degeneracy, 
however this does not yet explain the role of the j quan- 
tum number, and the chiral partners in the spectrum 
[lit . Another approach consists in considering a string 
with static ends, applicable to heavy quarks. In this case 
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the static quark-antiquark potential is extended from the 
hnear potential to, 



(3) 



where the pre-factor 2 of the radial excitation quantum 
number n is now the double of the pre-factor 1 in eq. ^ 
obtained with transverse modes only, f is a constant. 
This potential also fits correctly the lattice QCD exci- 
tations of the static quark-antiquark potential measured 
by Juge, Kuti and Morningstar [l^. 

We now briefly review how the potential of eq. ([3]) has 
been obtained by Buisseret, Mathieu and Semay, utilizing 
a simple 0, [21| constituent quark, antiquark and gluon 
hamiltonian, 



i=q,q,g j=q,q 
and applying einbeins, or auxiliary fields. 
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Einbeins were app lied to strings by Morgunov, Nefediev 
and Simonov [22l | to include the angular momentum of 
rigid strings in the quark confining potentials. Then 
Buisseret, Mathieu and Semay utilized the einbeins to 
also include the string excitations. The einbeins trans- 
form a more difficult problem into a harmonic oscillator- 
like Schodinger equation. Two independent oscillating 
modes can be separated, one in the usual r = — 
relative coordinate and another in the y = """^ '' — Vg 
orthogonal coordinate. This results in a spectrum, be- 
fore the einbein minimization, of, 
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(2ng + lg + 3/2) + — +H + Tr^ + ^+iy (6) 
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In the heavy quark case, similar to the static limit, 
/i ~ m » iig and the minimization of eq. ^ leads to 
the string excitations with static quarks as in eq. ([3]) . An 
important technical detail to arrive at eq. ^ is the pre- 
scription of Buisseret, Mathieu and Semay Sa 2tt(j to 
correct the einbein enhancement of the radial and angular 
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20l 21 1, of excited mesons with a linear po- 
tential. Solving the relativistic Schrodinger equation, the 
spectrum is led by M ^ ^/A-Kan + Sal [^. Morgunov, 
Nefediev and Simonov |22| showed that adding to the an- 
gular momentum of relativistic quarks the angular mo- 
mentum of a rigid and relativistic rotating string, the 
angular momentum contribution is corrected to M ^ 



y^2na{2n + I). On the other hand the einbein minimiza- 
tion leads to M ^y8a{2n + I). Therefore using the 
einbein, minimization together with the Scr -~^ 2'Ka pre- 
scription leads to the spectrum M — > y^27ra(2n+7) of 
the solution of the relativistic Schrodinger equation for 
a quark, an antiquark, and a non-vibrating string. The 
problem of the pre-factor 2 of the radial excitation re- 
mains, nevertheless that is not a result of the harmonic 
oscillator appearing with the einbeins, but a result of the 
linear potential. The difference between the spectra pro- 
duced with the harmonic oscillator and the linear poten- 
tials is the exponent of the principal quantum number 
in the spectrum. In particular a linear potential con- 
tributes to the spectrum as — > 2Tra{2n + I) to linear 
Regge spectra. This differs from the Harmonic Oscillator 
spectra M uj{2n + I). Also notice that PB Q tested 
enlarging the anzatse to the class of different power law 
confining potentials but this did not solve the problem of 
the principal quantum number, and it also may produce 
non-linear Regge trajectories. So one would better main- 
tain the potential linear, and look for another solution to 
the pre-factor 2 problem. 

Importantly, eqs. ([2]), Q and and the Regge be- 
haviour M ^2Tia{2n + I) of relativistic (but spinless) 
quarks suggest that a similar quantitatively correct angu- 
lar Regge slope a = {2iTa)~^ is possible with both quark 
and gluonic excitations. However these eqs. still lack the 
quark spin degrees of freedom, in particular they miss 
the insensitivity to chiral symmetry breaking observed 
in the large degeneracy of the excited meson spectrum. 
To solve this j quantum number problem, it is unavoid- 
able to work with confined and relativistic quarks with 
spontaneous chiral symmetry breaking. Notice that the 
spontaneous chiral symmetry breaking does not affect the 
Regge slopes, but it does correct the intercept, in partic- 
ular the TT meson is massless in the chiral limit, and the 
p meson mass is reproduced. The p sits in the lead- 
ing Regge trajectory and determines its intercept. High 
in the spectrum the good angular quantum number is 
not the orbital one / but the total one j. The neces- 
sary techniques were developed by Le Yaouanc, Oliver, 
Ono, Pene and Raynal [llj, and by PB and Ribeiro [l^ . 
and extended to hybrids by Llanes-Estrada, General and 
Cotanch [2^, . But in this letter we neglect the spin 
contribution of the gluon to the hamiltonian, since the 
fine structure of the static potential, also measured in 
Lattice QCD [2^, suggests that the spin Sg degree of 
freedom of the gluon only contributes at small distances, 
and Sg is not a leading term for excited mesons. The 
hamiltonian for the quarks is. 



H = 



■0t(.T) (mo/3 — id • V) V'(^) + 2-9^ 



— A*" 
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where it is standard in the Coulomb Gauge to use 
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a density-density colour octet potential V^^(r) = 
<5^iOi5i'0<5''''(— 3/16)F(r), with a funnel potential for the 
spatial function V{r) ~ ^ + f o + cr. Notice that for 
excited states we can neglect the Coulomb part j of the 
potential [2^ since it would have little effect in the ex- 
cited spectrum because the quark wavefunctions have a 
greater radius mean square. Moreover, when chiral sym- 
metry breaking occurs, the constant term vq is cancelled 
in the spectrum by an opposite term present in the quark 
self-energy. Thus for the study of quark degrees of free- 
dom in excited mesons it is sufficient to consider the lin- 
ear potential ar only. 

Here we go beyond these previous studies of mesons 
with chiral invariant quarks, including in this framework 
the gluon or string excitations, as in eq. (|4]), considering 
a vanishing current quark mass niq = 0. We also ap- 
ply the einbein technique to simplify the resulting three- 
body problem as in eq. (O. This allows us to sepa- 
rate the two coordinates r and y. In what concerns the 
relative quark-antiquark r coordinate, we address chiral 
symmetry breaking the bethe Salpeter equation for the 
quark and the antiquark with the framework already de- 
veloped for mesons. We now review the technical steps 
of this framework, before including the gluon or string 
degree of freedom, corresponding to the coordinate y. 

The boundstate equations are derived translating the 
relativistic invariant Dirac-Feynman propagators 111, in 
the quark and antiquark Bethe-Goldstone propagators 
0j l^rt - used in the formalism of non- relativistic quark 
models, 
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vs{0), (8) 



where it is convenient to define different functions of 
the constituent quark mass mc{k), (p{k) = arctan ^^IeIM ^ 
S{k) = sin(p(fc), C = cos(p(fc) and g{k) = 1 - S{k). 
In the non condensed vacuum, mc{k) equals the bare 
quark mass toq . In the physical vacuum, the constituent 
quark mass rndk) is a variational function which is de- 
termined by the mass gap equation. But here we do not 
detail it since the constituent mass is only crucial for 
the groundstate mesons, but not for the higher states in 
the spectrum where it can be neglected when compared 
to the momentum of the quark. The Salpeter-RPA equa- 
tions for a meson (a colour singlet quark-antiquark bound 
state) can be derived from the Lippman-Schwinger equa- 
tions for a quark and an antiquark, or replacing the prop- 
agator of eq. ([HI) in the Bethe-Salpeter equation. For a 
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[-2v'qV'i + W^lC^] [(S, ■ k){Sq ■ k) - • S,-] 



TABLE L The positive and negative energy spin-independent, 
spin-spin, spin-orbit and tensor potentials, computed exactly 
in the framework of the simple density-density harmonic os- 
cillator. For a detailed derivation, see reference 0, [23| . 



detailed derivation of the boundstate equations see ref- 
erence H^]. One gets. 



R+{k,P) 
R-\k,P) 

x{k,p) 



u\ki)x{k,P)v{k2) 

-M{P)-E{ki)-E(k2) 

vHk,)x{k,P)u{k2) 
-M{P)-E{k^)^E{k2) 
d^k' 

^V{k- k')[u{k[)R+{k\P)v\k',) 
+vik[)R-\k\P)uHk',)\ (9) 



where V{k) is the Fourier transform of the potential in 
eq. ([7]), and where /ci = fc + -j , A:2 = fc — -j and P 
is the total momentum of the meson. Eq. ^ includes 
both the Salpeter-RPA equations of PB et al. and 
of Llanes-Estrada et al. |28ll and the Salpeter equations 
of Le Yaouanc et al. [ll|]. Substituting x we get the 
equation for the positive energy i?"*" and negative energy 
i?^ radial wavefunctions. This results in four potentials 
V"^ with a = ±, respectively coupling = k R" to , 
in the boundstate Salpeter equation. 
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p- = Mp^ 
) p- = -Mp- 



(10) 

We now apply some educated approximations to arrive 
at an analytical expression for the excited meson spec- 
trum. As in eq. ([5]), we use einbeins and this avoids 
the problem of computing the matrix elements of the 
linear potential, actually we only need to know the ma- 
trix elements of the harmonic oscillator potential V^q 
as shown [3 in Table H and to minimize each energy in 
the spectrum with regards to p, v and pg. Importantly, 
the potentials V^'^ = V~~ and = V ^ include 

the usual spin-tensor potentials [1, [i^l , produced by the 
Pauli a matrices in the spinors of eq. ([8]) . They are de- 
tailed explicitly in Table |T1 Moreover we are interested 
in highly excited states, where both (r) and (k) are large 
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due to the virial theorem, and thus we consider the hmit 
where ^ as in PB, Cardoso, Van Cauteren and 
Llanes [29|. This imphes that in the potentials Hsted 
in Table [H we may simplify ip'{k) — *■ 0, C{k) — > 1 and 
G{k) = 1 — S{k) — > 1. We also notice that the relativistic 
equal time equations have twice as many coupled equa- 
tions as does the Schrodinger equation. But the negative 
energy component p~ is smaller than the positive energy 
component by a factor of the order of ^Ja jM. Thus 
when the meson mass M is large, and this is the case 
for the excited mesons, the negative energy components 
can be neglected and the Salpeter equation simplifies to a 
Schrodinger equation, where the relevant potential 
becomes. 
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the dependence on \? is traded by a dependence on = 
L2 + 2L-S + S^ The other remaining potentials are the 
spin- spin and tensor potentials, both much smaller than 
fc^ or than J^, and they simply contribute to a Coulomb- 
like potential p-. 

Including the gluonic coordinate y and Fourier trans- 
forming it as well, removing the Coulomb terms un- 
affecting the excited spectrum, we get simple Harmonic 
Oscillator terms. Substituting the Harmonic Oscillator 
matrix elements, we arrive at the meson spectrum. 
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before minimization, similar to eq. ^ except that the 
quark orbital angular momentum Iqq is replaced by the 
total quark angular momentum jqq. Eq p2p is hard to 
minimize analytically, except in two different limits, the 
limit of a large effective gluon mass, and the limit of a 
small number of gluon excitations. Then we get. 
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(13) 



where we denote Afqq ~ 2nqq+ jqq +3/2 and Mg = 2ng + 
Ig + 3/2. In the case of a vanishing effective gluon mass 
or of a significant Afg we have to resort to a numerical 
minimization of eq. (new einbein spectrum) (again with 
the 8cr 21X0 prescription) and we get the fit 



M 



^2TTa{Nqq+Ug) 



1.29 + 0.18:^^i— ^ 

Man + K 



Mqq + Afg 
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in terms of a dominant principal number Afqq + Afg , and 
a small contribution from a second number Afqq — Afg. 
Notice that the principal quantum number is in general 
much larger than the second number, due to the zero 
point energy of the Harmonic Oscillator. The simplest of 
all these cases is the one of a massless gluon and of few 
gluonic excitations, producing 



M 



^2tT(T {Afqq + Afg) 



(15) 



also leading to linear Reggc trajectories. 

Thus in any of the cases considered we find, as a 
good approximation that a principal quantum num- 
ber Afqq + Afg dominates the excited light meson spec- 
trum. We now show that the principal quantum number 
Afqq +Afg = 2nqq+jqq + 2ng + lg + 6 cau be simplified to 
a form n + j. This is quite interesting since it allows us 
to understand the experimental spectrum with both the 
quark quantum numbers and the gluon quantum num- 
bers! Obviously, the trajectories in Jtotai = Jgg+Lg have 
the correct slope of {2tt<t)^^. Thus we only need to show 
that the radial-like trajectories, i. e. the ones that main- 
tain jtotai and the parity P, also have the same slope, 
even in the less favourable case where the gluonic radial 
excitation pre- factor is 2 and not 1 as in eq. ([2]) . As in the 
radial trajectories of Bugg [l| , let us consider a radial tra- 
jectory starting with mesons which already have angular 
quark excitations jqq > 0. Due to the chiral symmetry in- 
sensitivity we start with a degenerate parity doublet, in- 
cluding a meson with jtotai — jgq+jg and P = + and an- 
other with jtotai = jqq+jg and P = -, i. e. respectively 
with j^ = and Jj~jjj;- The groundstate, in the lead- 
ing trajectory, has with no hybrid excitation, i. e. it has 
rig ~ and with Ig = Then, increasing Ig by one unit 
and summing it to the quarks angular momentum jqq 
with the triangular relation \jqq — lg\ < jtotai < jqq + Ig 
it is clearly possible to maintain the initial jtotai = jqq- 
Although the gluon angular momentum does invert the 
parity of the initial state, as depicted in Fig. [21 since we 
start with a parity doublet, this gluonic excitation 



Aj = , AM 



/r2 



2TTa 



(16) 



leads to another parity doublet. This doublet is a radial- 
like excitation of the one we started from. But, impor- 
tantly, the splitting for this radial-like excitation is 
only 27rcr. The next excitation, with AAf^ = 2 x 27rcr 
is easier to get since it can be intrinsically radial, with 
A{2nqq + cug) = 2. And so forth, with the interplay 
of Jgg and Lg we are able to produce new radial-like 
excitations, in the middle of the purely radial excita- 
tions. Thus the Regge slope for the radial-like excitations 
a = {2TTa)~^ is the same one of the angular excitations. 

To conclude, we address a novel large degeneracy ob- 
served in quantum physics, with a principal quantum 
number n-\-j in the excited light meson spectrum. With 
the interplay of quark excitations insensitive to chiral 
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j =jqa+jg-2, or 2nqg^+ng=2 
j=jaa+lg' lg=o • nqa+ng=0 



J 



J 



FIG. 2: Mass squared M^, as a function of , of radial-like 
excitations in 2na units starting from the groundstates in the 
total angular momentum j shell. 



symmetry and of gluon string excitations we show that 
the degeneracy may be explained. We simplify the quark 
with strings problem to a 3-body problem with two chiral 
invariant quarks and a gluon. An exact numerical solu- 
tion of this three-body problem would require a similar 
effort to the one necessary to study other three-body sys- 
tems of chiral invariant quarks and gluons [1^, [13] , e. g. 
9 dimensional spatial integrals, the sum of 6 spin indices, 
and the use of a variational basis with hundreds of multi- 
dimensional orthogonal functions, implying man-years of 
coding and debugging together with the use of a super- 
computer for several months. Thus for this first study, 
and to to arrive at simple and clear analytical mass for- 
mulas, we use the einbein technique. As a by-product of 
our research, we find that the einbein technique is appli- 
cable to the quark models with chiral symmetry breaking. 
For instance, eq. (|15p . with the quark degrees of freedom 
only, should provide a good approximation to the exact 
numerical solution of Llanes-Estrada and Cotantch [2^ 
and of Glozman and Wagenbrunn 0]. 

Note added in proof - Very recently Bugg re- 
analysed Cristal Barrel data at LEAR (sij and his spec- 
trum shows clearly the large degeneracy, where the angu- 
lar excitations have the same scale as the radial-like ex- 
citations. However he found no evidence for parity dou- 
blets in the leading Regge Trajectory, but only in daugh- 
ter trajectories. More experimental, computational and 
theoretical works on the light meson spectrum are urgent. 

PB tanks the colleagues at the 2009 Schladming Win- 
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plete this paper. PB was supported by the FCT 
grants POCI/FP/81913/2007, POCI/FP/81933/2007 
and CERN/FP/83582/2008. 
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